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ABSTRACT 

Context. Massive young stellar objects (MYSOs) with hot cores are classic sources of 
complex organic molecules. The origins of these molecules in such sources, as well as the 
small- and large-scale differentiation between nitrogen- and oxygen-bearing complex 
species, are poorly understood. 

Aims. We aim to use complex molecule abundances toward a chemically less explored 
class of MYSOs with weak hot organic emission lines to constrain the impact of hot 
molecular cores and initial ice conditions on the chemical composition toward MYSOs. 
Methods. We use the IRAM 30m and the Submillimeter Array to search for complex or¬ 
ganic molecules over 8-16 GHz in the 1 mm atmospheric window toward three MYSOs 
with known ice abundances, but without luminous molecular hot cores. 

Results. Complex molecules are detected toward all three sources at comparable abun¬ 
dances with respect to CH 3 OH to classical hot core sources. The relative importance of 
CH 3 CHO, CH 3 CCH, CH 3 OCH 3 , CH 3 CN, and HNCO differ between the organic-poor 
MYSOs and hot cores, however. Furthermore, the N-bearing molecules are generally 
concentrated toward the source centers, while most O- and C-bearing molecules are 
present both in the center and in the colder envelope. Gas-phase HNCO/CH 3 OH ratios 
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are tentatively correlated with the ratios of NH 3 ice over CH 3 OH ice in the same lines 
of sight, which is consistent with new gas-grain model predictions. 

Conclusions. Hot cores are not required to form complex organic molecules, and source 
temperature and initial ice composition both seem to affect complex organic distribu¬ 
tions toward MYSOs. To quantify the relative impact of temperature and initial condi¬ 
tions requires, however, a larger spatially resolved survey of MYSOs with ice detections. 

Key words. Astrochemistry, ISM : molecules, ISM : abundances 


1. Introduction 

Organic molecules containing more than six atoms, the so-called complex organics (Herbst 
& van Dishoeck 2009), are commonly found in the warm and dense gas (T > 100 K, n > 
10 6 cm~ 3 ) around young stellar objects (YSOs), so-called molecular hot cores (e.g., Blake 
et al. 1987; Cazaux et al. 2003; Fuente et al. 2005). Abundances and abundance ratios of 
complex organics are found to vary substantially between (Helmich & van Dishoeck 1997) 
and within YSOs (e.g., Wyrowski et al. 1999). This suggests that formation and destruc¬ 
tion routes are highly environment specific and that there is a sensitive dependence of 
the complex organic chemistry on chemical and physical initial conditions. In addition, 
different filling factors of the warm gas should play a role if there the complex organic 
products of cold and hot chemistry differ there. 

The potential environmental dependencies and chemical memories lead to complex 
organics having a great potential as probes of the current and past physical and chemical 
conditions where they are found (Nomura & Millar 2004). Their potential utility is further 
increased by the fact that most complex organic molecules present large numbers of lines, 
spanning most excitation conditions found in space. Complex molecules are also of high 
interest for origins of life theories since they are the precursors of even more complex 
prebiotic material (Ehrenfreund & Charnley 2000). Using molecules as probes of physical 
conditions and advancements in prebiotic evolution from organics both rely on a detailed 
understanding of complex organic chemistry. The formation and destruction mechanisms 
and rates of most complex organics are, however, poorly constrained. 

The formation of organic molecules around massive YSOs (MYSOs) was first thought 
to proceed through gas phase reactions in dense hot cores, following evaporation of ice 
grain mantles (e.g., Charnley et al. 1992). Recent laboratory experiments and modeling 
efforts point now toward a more complicated sequential scenario that relies to a greater 
extent on surface formation routes on submicron-sized dust particles. Herbst & van 
Dishoeck (2009) classify complex organic molecules in terms of generations according to 
the following scenario. In interstellar clouds and in the deeply embedded early phases of 
star formation, atoms and molecules accrete or form on the surface of dust grains, building 

Article number, page 2 of 28 



Edith C. Fayolle et al.: New constraints on the origins of organics 


up an icy mantle of simple species like H 9 O, CH 4 , and NH 3 (Tielens & Hagen 1982). This 
icy mantle is processed at low temperature by atoms, which can diffuse even at the low 
temperatures in cloud cores, creating the zeroth generation of organic molecules. A good 
example of these species is CH 3 OH, which is efficiently formed at low temperature by 
the hydrogenation of CO ice (Watanabe & Kouchi 2002; Watanabe et al. 2003, 2004; Fuchs 
et al. 2009; Cuppen et al. 2009). First-generation complex organics form when heating the 
cold envelope up by the increasing luminosity of a central YSO and is due to a combina¬ 
tion of photoprocessing of the ice resulting in radical production and a warming up (20 
to 100 K) of the grains, thereby enhancing the mobility of radicals and molecules (Gar- 
rod et al. 2008; Oberg et al. 2010). When the icy grains move inward and reach a region 
warmer than 100 K, the icy mantle evaporates, bringing the zeroth- and first-generation 
organics into the gas phase, where additional chemical reactions give rise to the formation 
of the second-generation complex organics (e.g., Charnley et al. 1992; Doty et al. 2002; Viti 
et al. 2004). 

In the proposed scenario of complex molecule formation, the initial ice mantle plays 
a critical role. The exact composition of this ice may therefore have a strong effect both 
on the product composition of formed organics and on their overall formation efficiency. 
Garrod et al. (2008) and Oberg et al. (2009) find, for example, that CH 3 OH ice is a key 
starting point for most complex organic formation. Rodgers & Charnley (2001) used a hot 
core chemistry model to show that the relative amount of NH 3 in the ice has a large impact 
on the CH 3 CN/CH 3 OH protostellar abundance ratio. Observationally testing these rela¬ 
tionships would provide key constraints on the formation pathways of complex organic 
molecules. 

Isolated MYSOs with warm inner envelopes are good laboratories for testing this hy¬ 
pothesis as these sources are bright enough to observe a wide variety of organics and 
some of them present ice features from the cold outer protostellar envelope (Gibb et al. 
2004). Sources presenting both complex gas and ice features are, however, rare as the 
sources need to be evolved to possibly display a bright hot core chemistry accessible to 
current observational facilities and young enough such that the ice material has not been 
completely consumed by accretion, warm up, and envelope dispersal. In the massive 
YSO sample studied by Bisschop et al. (2007), only three hot cores present ice spectra (see 
Table 1). Such a small number prevents any analysis of the correlation between ice and 
gas content and justifies our search for other objects that display both ice features and gas 
phase organics. 

To extend the sample of sources with both complex organics and ice observations, we 
look for gas phase organics species around non-hot core MYSOs (absence or low-level of 
hot CH 3 OH emission) that also have ice observations available from the literature. These 
sources are called from now on organic-poor MYSOs (poor in lines of organic molecules). 
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Complex molecule observations in such objects may additionally shed light on the con¬ 
ditions under which different kinds of complex molecules can form, i.e. which molecules 
require the presence of a hot core to be abundant. 


Massive objects NGC7538 IRS9, W3 IRS5 , and AFGL490 have been observed in the 
mid-infrared by the Infrared Space Observatory (ISO) and analyzed systematically for ice 
abundances by Gibb et al. (2004) and references therin. NGC7538 IRS9 is a 6 x 10 4 Lq lu¬ 
minous object located in Perseus. It is close to hot core source NGC7538IRS1 and displays 
at least three bipolar outflows, evidence for accretion (Sandell et al. 2005), and a hot com¬ 
ponent close to the central object. W3 IRS5 is associated with five YSOs, two of which are 
massive (van der Tak et al. 2005; Megeath et al. 2005; Rodon et al. 2008; Chavarria et al. 
2010). It has a luminosity of 17 x 10 4 Lq and presents strong S-bearing molecular lines 
(Helmich et al. 1994). AFGL490 is a very young medium-mass YSO of 4.6 x 10 3 Lq, in 
transition to a Flerbig Be star, which drives a high-velocity outflow (Mitchell et al. 1995) 
and shows evidence of a rotating disk (Schreyer et al. 2006). Since these sources are con¬ 
sidered to potentally be at an earlier evolutionary stage than typical hot-core sources, it is 
difficult to predict the chemical complexity and the spatial emission of the organics that 
could be observed in these sources. 


In this study we use a combination of single-dish IRAM 30m data and spatially re¬ 
solved observations from the Submillimeter Array 1 (SMA) to search for organic molecules 
around these three MYSOs and report on their complex organic abundances in the cool 
protostellar envelope and in a warmer region closer to the star. A subset of these data 
was used in the Oberg et al. (2013) to study the detailed radial distribution of molecules 
in NGC7538 IRS9, while the present study focuses on the overall detection rate of organ¬ 
ics in these organic-poor sources, and on how they compare with ice abundances and 
traditional hot core chemistry. The paper is organized as follows. The observations are 
described in Section 2, and the results of the line analysis are shown in Sections 3.1, 3.2, 
and 3.3. The chemistry in our sample is compared to the chemistry in traditional hot-core 
sources in Section 3.4. Section 3.5 presents correlation studies between ice and gas column 
densities and abundances, testing the impact of initial ice compositions on the complex 
chemistry. A discussion of the use of these line-poor sources to underpin the origins of 
complex chemistry is presented in Section 4, which is followed by the conclusions of this 
study. 
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Table 1: Source characteristics and ice abundances. The sources observed in this study are 
in boldface, the others are from Bisschop et al. (2007). 


Source 

a:( 2000 ) 

<5(2000) 

d 

kpc 

L 

10 4 L o 

Nh ? o 

10 r/ cm -2 

CH 3 OH 

X [%] (/N h ,o ) 
ch 4 nh 3 

OCN~ 

NGC7538 IRS9 

23:14:01.6 

+61:27:20.4 

2.7 

3.5 

70 

4.3 + 0.6 

2 + 0.4 

15 + 2.7 

1.7 + 0.5 

W3 IRS5 

02:25:40.5 

+62:05:51.3 

2.0 

17 

51 

<3.3 

<1.3 

<5.7 

<0.23 

AFGL490 

03:27:38.7 

+58:47:01.1 

1.4 

0.46 

6.2 

11 + 4 

<2.4 

<16 

< 1.2 

W33A 

18:14:38.9 

-17:52:04.0 

3.8 

5.3 

110 

15 + 5 

1.5 + 0.2 

15 + 4 

6.3 + 1.9 

AFGL2591 

20:29:24.6 

+40:11:19.0 

3.3 

18 

12 

14 + 2 

<2.7 

<2.3 

- 

NGC7538 IRS1 

23:13:45.4 

+61:28:12.0 

2.4 

15 

22 

<4 

1.5 + 0.5 

<17 

<0.5 

Orion IRc2 

05:35:14.3 

-05:22:31.6 

0.4 

1.0 

24.5 

10 + 3 

- 

- 

2 + 0.6 

G24.78 

18:36:12.6 

-07:12:11.0 

7.7 

1.2 

- 

- 

- 

- 

- 

G75.78 

20:21:44.1 

+37:26:40.0 

1.9 

19 

- 

- 

- 

- 

- 

NGC6334IRS1 

17:20:53.0 

-35:47:02.0 

1.7 

11 

- 

- 

- 

- 

- 


2. Observations and analysis 

2.1. Observations 

The MYSOs NGC7538 IRS9, W3 IRS5 , and AFGL490 located in Perseus NGC7538 at 2.7 
kpc, in Perseus W3 at 2.0 kpc, and in Camelopardalis OBI at 1.4 kpc respectively (see Table 
1) were observed with the IRAM 30m and the Sub millimeter Array. The three sources 
were observed with the IRAM 30m telescope on February 19-20, 2012 using the EMIR 
230 GFIz receiver and the new FTS backend. At these frequencies the IRAM 30m beam is 
~10 // . The two sidebands cover 223-231 GFIz and 239-247 GFIz at a spectral resolution of 
~0.2 km s” 1 and with a sideband rejection of -15dB (Carter, M. et al. 2012). We checked 
the pointing every one to two hours and found to be accurate within 7!' to 3". 

Focus was checked every four hours and generally remained stable through most of 
the observations; i.e., corrections in the range of 0.2-0.4 were common, but a correction 
of 0.7 was required once. We acquired spectra in both position-switching and wobbler¬ 
switching modes. The resulting spectra had similar relative line intensities, indicative of 
no emission in the wobbler-off position. The wobbler-switching mode was considerably 
more stable, and we used these data alone for the quantitative analysis. The weather 
during the observations was excellent and the T 225 GHZ varied between 0.05 and 0.15. 
We converted the raw IRAM spectra to main beam temperatures and fluxes using for¬ 
ward and beam efficiencies and antenna temperature to flux conversion values listed at 
www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies. The spectra were reduced us¬ 
ing CLASS 2 . A linear baseline was fitted to each 4 GFIz spectral chunk using four to seven 
windows. The individual scans were baseline-subtracted and averaged 3 . The absolute 
flux scale of the lines were then set using calibrated SMA data as outlined in detail by 
Oberg et al. (2013). 

1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory 
and the Academia Sinica Institue of Astronomy and Astrophysics. It is funded by the Smithsonian 
Institute and the Academia Sinica. 

2 CLASS website: http://www.iram.fr/IRAMFR/GILDAS 

3 Reduced data are available through the dataverse network at http: //dx. doi. org/10.7910/DVN/ 
26562 
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Fig. 1: Image of the CH 3 CN emission using the 13o - 12q at 239.138 GHz line acquired by 
the SMA for the massive young stellar objects NGC7538 IRS9, W3IRS5, and AFGL490 tar¬ 
geted in this study. The black contour presents the 50 % line intensity, and the synthesized 
beam is shown in white at the bottom left. A 2" radius mask used to extract the spectra is 
overplotted in dashed red line. Images for the hot core source NGC7538IRS1 is presented 
as well. The latter source has been through the same program as the three other sources. 


SMA observations were acquired in the compact and extended array configurations. 
The data in the compact configuration were taken on 15 October 2011 for all sources and 
with seven antennas, resulting in baselines between 16 m and 77 m. The data in the 
extended configuration were obtained using eight antennas, resulting in 44 m to 226 m 
baselines and were acquired on 29 July 2011 for W3 IRS5 and AFGL490 and on the 15th 
of August 2011 for NGC7538 IRS9. We set-up the SMA correlator to obtain a spectral res¬ 
olution of ~1 km s. 1 using 128 channels for each of the 46 chunks covering 227-231 GHz 

in the lower sideband and 239-243 GHz in the upper sideband. The T 225 GHz was 0-09 on 
29 July, 0.1 on 15 August, and 0.07 on 15 October 2011 4 . 

We used the MIR package 5 to perform the first data reduction steps (flux calibration 
and continuum subtraction). Absolute flux calibration is done with Callisto. The band¬ 
pass calibrators 1924-292 and 3c84 were used for the compact observations, and 3c454.3 
and 3c279 were used to calibrate 29 July and 15 July observations, respectively. The 
quasars 0014+612 and 0102+584 were used as gain calibrators for NGC7538 IRS9, and 
0244+624, 0359+509, and we used 0102+584 for W3 IRS5 and AFGL490. The compact 
and extended data were combined for each source with MIRIAD 5 using natural or robust 
weighting, depending on the data quality, which resulted in synthesized beam sizes of 
2.0" x 1.7" for NGC 7538 IRS 9, 2.2" x 2.8" for W3 IRS5, and 2.3" x 2.9" for AFGL490. 

2.2. Spectral extraction and rms 

Both the IRAM and SMA data were frequency-calibrated using the bright 5-4 CH 3 OH 
ladder around 241.7 GHz, correcting for the intrinsic velocity of the different sources. We 
extracted the SMA spectra using a 2 /, -radius mask around the continuum phase center of 
each source. The mask dimension was chosen to encompass a majority of the CH 3 CN line 

4 Observations are available on the SMA archive website herhttp://www.cfa.harvard.edu/cgi- 
bin/sma/smaarch.pl 

5 MIR website: http://www.cfa.harvard.edu/~cqi/mircook.html 
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emission at 239.318 GHz that can be associated with a core component, as shown in Figure 
1. We selected a 2" mask size based on a combination of theory and data inspection, i.e. 
the optimal mask size should include all the hot emission and exclude as much as possible 
of the cold envelope emission. In all sources, the selected mask size should be bigger than 
the 100 K radius and thus incorporate all emission associated with a potential hot core. 
Toward NGC 7538 IRS9 and AFGL 490, where the 100 K radius should be smaller than 2" 
, smaller masks were also explored to more exclusively trace the T > 100 K region, but the 
resulting spectra had generally too low signal-to-noise ratio to be useful for a quantitative 
analysis. Some colder chemistry contribution to the SMA spectra in these sources cannot, 
thus, be excluded a priori, but the succeeding analysis (see below) demonstrated that the 
emission is indeed dominated by hot gas. 

The rms for the IRAM and SMA observations of each source was derived in a line 
free region of several hundred channels: the 229.37-229.445 GHz region for the lower 
side band and the 240.7 - 240.75 GHz region for the upper side band. The rms derived 
for the IRAM observations is between 15 and 20 mK, which is lower than any previous 
millimeter observations for these sources. For the SMA data, the rms for the lower side 
band is ~70 mK, and~100 mK for the upper side band. 

3. Results 

3.1. Line identification and characterization 

Figure 2 shows the IRAM 30m 239-243 GHz spectra for the three targeted line-poor MYSOs 
and the hot-core source NGC7538 IRS1. The organic-poor MYSOs have, as expected, a 
lower line density, but also many line coincidences with the hot core. Of the lines from 
complex organics listed by Bisschop et al. (2007) found in their sample hot-core sources, 
CH 3 OH, CH 3 CN, CH 3 CCH, HNCO, CH 3 OCH 3 , and CH 3 CHO lines were identified in at 
least one of the organic-poor MYSOs using the splatalogue catalog tool 6 and the CDMS 7 
and the JPL 8 spectral databases (Muller et al. 2001; Pickett et al. 1998). All available lines 
in the observed spectral range were used for the quantitative analysis except for CH 3 OH 
where we only used the lines from the 5-4 ladder to simplify the excitation analysis. 

We fitted the identified lines with a Gaussian function in IDL using the routine 'gauss- 
fit' for isolated lines and 'mpfitfun' when a multiple Gaussian fit was required because 
of overlapping lines. A local baseline component was added to the fits when needed, 
and the presented uncertainties were output by the fitting routines. We calculated 3 cr up¬ 
per limits using an average FWHM for the different sources. Unresolved multiplets were 
treated in one out three ways depending on the nature of the overlapping lines: 1) If one 

6 Splatalogue website: http://www.cv.nrao.edu/php/splat/ 

7 CDMS website: http://www.astro.uni-koeln.de/cdms 

8 JPL database website: http://spec.jpl.nasa.gov/ 
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Fig. 2: 239-243 GHz spectral window from the IRAM 30m displaying emission lines for 
typical hot core source NGC 7538 IRS land weak line MYSOs NGC7538 IRS9, W3IRS5, 
AFGL490. The star-marked lines are CO ghost lines consistent with the sideband rejection 
for each source. 

of the possible contributing lines had a very low Einstein coefficient or high upper en¬ 
ergy level and/or is not likely to be detected based on non-detections of the same species 
in other frequency ranges, then it was assumed to not contribute significantly and was 
not included in the fit. 2) If the lines came from the same species and the upper energy 
level and Einstein coefficients were identical or close to identical, then the degeneracies 
were added and the feature was treated as a single line, 3) if none of the two previous 
conditions were met, we did not include the multiplet in the analysis. 

Line upper energy levels, Einstein coefficients, degeneracies, and quantum numbers 
from the Splatalogue are listed with the derived line fluxes and FWHM in Table 2 for 
CH 3 OH from the single-dish observations, in Table 3 for CH 3 OH from the SMA spectra, 
in Table 4 for CH 3 CN from IRAM, in Table 5 for CH 3 CN from the SMA, in Table 6 for 
CH 3 CCH from IRAM, in Table 7 for CH 3 CCH from the SMA, and in Table 8 for HNCO, 
CH 3 OCH 3 , and CH 3 CHO. Only the lines with an Einstein coefficient logarithm higher 
than —4.5 and their upper level energy below 400 K are displayed in the tables. Due to 
the high line density for CH 3 OCH 3 and CH 3 CHO, only the lines with an upper energy 
level below 200K are shown for these species. No other complex molecules were detected 
toward any of the sources. For molecules with weak lines, we only used the IRAM data 
since the SMA observations have lower spectral resolution and signal-to-noise ratio. 

3.2. Spatial origin of the line emission 

Figures 3, 4, and 5 present the line fluxes of key molecules from both the single-dish and 
SMA observations toward the three MYSOs. The IRAM beam is 6.2 to 7.2 times larger 
than the SMA mask ((IRAM radius at 227-243 GHz: 5-5.4 // ) 2 /(SMA mask radius: 2 // ) 2 ). 
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Table 2: CH 3 OH lines data from IRAM 30m spectra. 


Freq 

(GHz) 

Eup 

(K) 

log A 

gu 

Transition 

NGC7538IRS9 
$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3 IRS5 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

AFGL490 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

239.746 

49.1 

-4.25 

11 

5i,s - 4 j, 4 A + 

20.7+2.3 

3.70+ 0.03 

6.1+0.8 

2.32+ 0.04 

4.6+0.7 

4.4+ 0.2 

241.700 

47.9 

-4.22 

11 

5o,5 -4 0< 4 E 

30.8+3.3 

3.47+ 0.02 

7.1+0.9 

2.43+ 0.04 

6.7+0.9 

3.53+ 0.07 

241.767 

40.4 

-4.24 

11 

5-1,5 -4-1,4 E 

60.7+6.3 

3.37+ 0.01 

11.2+1.3 

2.98+ 0.03 

16.7+1.9 

3.11+ 0.02 

241.791 

34.8 

-4.22 

11 

5o,5 -4(3,4 A 

69.2+7.1 

3.37+ 0.01 

12.6+1.5 

2.91+ 0.03 

20.3+2.2 

3.09+ 0.02 

241.807 

115.2 

-4.66 

22 

5 4 - 4 4 A- 

2.3+0.5 

4.5+ 0.4 

1.6+0.3 

1.58+ 0.09 

<1.4 

- 

241.813 

122.7 

-4.66 

11 

5-4,2 -4-4,1 E 

1.1+0.4 

3.6+ 0.6 

1.2+0.3 

1.7+ 0.2 

<1.4 

- 

241.830 

130.8 

-4.66 

11 

54,1 -44,o E 

<1.1 

- 

1.3+0.3 

2.0+ 0.2 

<1.4 

- 

241.833 

84.6 

-4.41 

22 

5 3 -4 3 A- 

13.3+1.6 

4.61+ 0.06 

5.5+0.7 

2.18+ 0.04 

3.6+0.6 

4.9+ 0.2 

241.842 

72.5 

-4.29 

11 

52,4 -42,3 A 

10.1+1.4 

5.5+ 0.1 

3.2+0.5 

1.90+ 0.08 

2.6+0.6 

6.8+ 0.5 

241.844 

82.5 

-4.41 

11 

53,2 - 43,1 E 

2.4+0.5 

3.3+ 0.3 

3.2+0.5 

2.2+ 0.1 

1.5+0.3 

3.50+ 0.02 

241.852 

97.5 

-4.41 

11 

5-3,3 -4-3,2 E 

3.2+0.6 

5.5+ 0.4 

2.2+0.4 

2.1+ 0.1 

1.5+0.4 

6.5+ 0.8 

241.879 

55.9 

-4.22 

11 

5i,4 -4i,3 E 

23.5+2.6 

3.64± 0.02 

6.7+0.9 

2.58+ 0.05 

5.2+0.7 

4.1+ 0.1 

241.888 

72.5 

-4.29 

11 

52,3 -42,2 A + 

8.5+1.1 

4.28+ 0.09 

3.8+0.5 

2.04+ 0.06 

2.5+0.5 

5.2+ 0.3 

241.904 

60.7 

-4.29 

11 

5-2,4 -4-2,3 E 

15.9+1.7 

3.66+ 0.02 

5.0+0.6 

2.57+ 0.03 

3.9+0.5 

3.85+ 0.07 

241.905 

57.1 

-4.30 

11 

52,3 -42,2 E 

15.9+1.7 

3.66+ 0.02 

5.0+0.6 

2.57+ 0.03 

3.9+0.5 

3.85+ 0.07 


Table 3: CH 3 OH lines extracted from SMA observations with a 2 /, -radius mask. 


Freq 

(GHz) 

E up 

(K) 

logA 

gif 

Transition 

NGC7538 IRS9 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3IRS5 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

AFGL490 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

239.746 

49.1 

-4.25 

11 

5i,s -4i, 4 A+ 

7.6+1.1 

4.6+ 0.2 

2.3+0.5 

1.9+ 0.2 

3.0+0.7 

4.6+ 0.5 

241.700 

47.9 

-4.22 

11 

5o,5 - 4o,4 E 

7.3+1.0 

4.1+ 0.2 

2.8+0.6 

2.4+ 0.2 

3.7+0.8 

5.9+ 0.5 

241.767 

40.4 

-4.24 

11 

5-1,5 - 4-1,4 E 

8.2+1.1 

3.5+ 0.1 

3.2+0.6 

2.3+ 0.2 

1.8+0.5 

2.6+ 0.3 

241.791 

34.8 

-4.22 

11 

5o,5 - 4(3,4 A 

7.1+1.0 

3.16+ 0.09 

3.5+0.7 

2.9+ 0.2 

2.4+0.6 

4.9+ 0.7 

241.807 

115.2 

-4.66 

22 

5 4 - 4 4 A- 

2.1+0.7 

6+ 1 

1.4+0.5 

2.4+ 0.5 

<1.1 

- 

241.813 

122.7 

-4.66 

11 

5-4,2 - 4-4,1 E 

1.0+0.4 

3.6+ 0.9 

1.0+0.5 

2.8+ 0.8 

<1.1 

- 

241.830 

130.8 

-4.66 

11 

54,1 - 4 4 ,o E 

<0.8 

- 

<1.1 

- 

<1.1 

- 

241.833 

84.6 

-4.41 

22 

5 3 -4 3 A- 

6.4+1.1 

4.7+ 0.3 

3.1+0.7 

2.9+ 0.3 

1.8+0.6 

4.6+ 0.8 

241.852 

97.5 

-4.41 

11 

5-3,3 - 4-3,2 E 

1.9+0.6 

3.4+ 0.5 

1.4+0.5 

1.5+ 0.4 

<1.1 

- 

241.879 

55.9 

-4.22 

11 

5l,4 - 4l,3 E 

5.5+0.9 

3.4+ 0.2 

2.6+0.6 

2.3+ 0.2 

1.9+0.6 

4.7+ 0.7 

241.888 

72.5 

-4.29 

11 

52,3 - 4 2,2 A+ 

4.6+0.8 

5.2+ 0.3 

1.6+0.5 

1.6+ 0.2 

1.8+0.7 

7+2 

241.904 

60.7 

-4.29 

11 

5-2,4 - 4-2,3 E 

5.3+0.8 

4.5+ 0.1 

1.7+0.3 

2.2+ 0.2 

1.7+0.4 

5.2+ 0.5 

241.905 

57.1 

-4.30 

11 

52,3 - 4 2,2 E 

5.3+0.8 

4.5+ 0.1 

1.7+0.3 

2.2+ 0.2 

1.7+0.4 

5.2+ 0.5 


Table 4: CH 3 CN lines data from IRAM 30m spectra. 


Freq 

(GHz) 

E up 

(K) 

logA 

g« 

Transition 

NGC7538IRS9 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3 IRS5 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

AFGL490 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

239.023 

258.9 

-3.00 

54 

13s -12 5 

1.6+0.5 

9+2 

<1.0 

- 

<1.5 

- 

239.064 

194.6 

-2.97 

54 

13 4 -12 4 

2.0+0.7 

9+2 

0.9+0.3 

6 + 2 

1.2+0.4 

6+2 

239.096 

144.6 

-2.95 

108 

13 3 - 12s 

5.3+0.8 

6.0+ 0.3 

2.2+0.4 

2.2+ 0.2 

2.4+0.6 

8 + 1 

239.120 

108.9 

-2.94 

54 

13 2 -12 2 

4.5+0.7 

5.7+ 0.3 

1.7+0.4 

2.8+ 0.3 

1.5+0.5 

6+2 

239.133 

87.5 

-2.93 

54 

13i - 12i 

5.8+0.8 

4.6+ 0.2 

2.4+0.4 

2.4+ 0.2 

2.1+0.5 

4.9+ 0.5 

239.138 

80.3 

-2.93 

54 

13 0 - 12q 

7.4+1.0 

5.2+ 0.2 

2.6+0.4 

2.1+ 0.2 

2.4+0.5 

5.4+ 0.5 


That most emission lines in these figures do not display a factor of six or seven difference 
between the IRAM 30m and SMA spectra demonstrates a non-uniform emission across 
the object. Some emission line fluxes, most notably CH 3 CN, are similar (within a factor 
of two) between the IRAM 30m and SMA spectra, indicating of a large contribution from 
unresolved emission at the source center. In contrast, little or no CH 3 CCH flux from the 
IRAM is recovered by the SMA, which indicates extended emission. The fact that some 
CH 3 CCH IRAM 30m fluxes are more than 6 .2-7.2 x higher than the corresponding SMA 
fluxes is explained by spatial filtering of large-scale emission and / or off-centered emis- 
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Table 5: CH 3 CN lines data from SMA spectra. 


Freq 

(GHz) 

E up 

(K) 

logA 

g« 

Transition 

NGC7538IRS9 
$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3 IRS5 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

AFGL490 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

239.064 

194.6 

-2.97 

54 

13 4 -12 4 

2.9+0.9 

7.6+ 1.2 

< 0.8 

- 

< 1.9 

- 

239.096 

144.6 

-2.95 

108 

13 3 - 12s 

5.1+1.0 

5.6+ 0.4 

2.9+0.9 

6 .8+ 1.2 

< 1.9 

- 

239.120 

108.9 

-2.94 

54 

13 2 -12 2 

4.4+1.1 

9.6+ 1.0 

1.2+0.5 

2.9+ 0.9 

1 .8+0.8 

6+2 

239.133 

87.5 

-2.93 

54 

13i - 12j 

5.6+1.3 

6.9+ 0.7 

3.1+0.6 

4.4+ 0.4 

1.8+0.9 

6+2 

239.138 

80.3 

-2.93 

54 

13 0 -12 0 

5.2+1.2 

6.3+ 0.6 

2.1+0.5 

2.7+ 0.4 

1.8+0.9 

6+2 


Table 6 : CH 3 CCH lines data from IRAM 30m spectra. 


Freq E up logA g„ Transition NGC7538IRS9 W3IRS5 AFGL490 


(GHz) 

(K) 




I FdV 

(Jykms J ) 

FHWM 
(km s -1 ) 

I FdV 

(Jy km s : ) 

FHWM 
(km s -1 ) 

I FdV 

(Jy km s - 1 ) 

FHWM 
(km s -1 ) 

239.088 

346.1 

-4.07 

16 

14 6 - 13 6 

<1.1 

- 

<0.7 

- 

<1.0 

- 

239.179 

201.7 

-4.88 

58 

14 4 - 13 4 

1 .0+0.3 

3.4+ 0.6 

0.9+0.3 

2.5+ 0.6 

<1.0 

- 

239.211 

151.1 

-4.00 

16 

14a - 13 3 

4.7+0.7 

3.1+ 0.2 

4.0+0.5 

2 .00+ 0.06 

1.5+0.3 

3.5+ 0.4 

239.234 

115.0 

-4.85 

58 

14 2 - 13 2 

5.1+0.7 

2.71+ 0.08 

3.6+0.5 

1.98+ 0.07 

1.5+0.4 

3.0+ 0.3 

239.248 

93.3 

-4.84 

58 

14i - 13, 

9.0+1.1 

2.92+ 0.05 

5.8+0.7 

2.2+ 0.05 

2.7+0.4 

2.3+ 0.1 

239.252 

86.1 

-4.84 

58 

14 0 -13 0 

10 .2+1.2 

2.74+ 0.04 

6.4+0.8 

2.12+ 0.04 

3.7+0.5 

2.9+ 0.1 


Table 7: CH 3 CCH lines data from SMA spectra. 


Freq 

(GHz) 

E up 

(K) 

logA 

§» 

Transition 

NGC7538IRS9 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3 IRS5 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

AFGL490 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

239.088 

346.1 

-4.07 

16 

14e -136 

<0.9 

- 

<0.7 

- 

<1.0 

- 

239.179 

201.7 

-4.88 

58 

14 4 -13 4 

<0.9 

- 

0.9+0.3 

2.5+ 0.6 

<1.0 

- 

239.211 

151.1 

-4.00 

16 

14 3 -13 3 

1 .8+0.6 

3.2+ 0.5 

<1.1 

- 

1 .0+0.6 

4.2+ 1.5 

239.234 

115.0 

-4.85 

58 

14 2 -13 2 

1.3+0.5 

3.5+0.7 

<1.1 

- 

<1.0 

- 

239.248 

93.3 

-4.84 

58 

14j -13 4 

2.3+0.7 

3.5+ 0.5 

<1.1 

- 

<1.0 

- 

239.252 

86.1 

-4.84 

58 

14 0 - 13q 

2.4+0.7 

2.9+ 0.4 

<1.1 

- 

1 .6+0.8 

5.9+ 1.9 


sion. CH 3 OH lines display a mixed behavior: lines with higher upper energies show 
more overlap between the IRAM 30m and SMA spectra than the colder lines. The IRAM 
30m and SMA line fluxes for the llo, 11 - 10o,io HNCO line at 241.774 GHz are close for 
NGC7538IRS9 and similar for W3IRS5, but none of the IRAM 30m flux is recovered by the 
AFGL490 SMA observations. Based on these lines, HNCO emission appears to be coming 
from both the core and the envelope of NGC7538 IRS9, from the core of W3 IRS5 alone, 
and from the envelope of AFGL490. This source-to-source difference could partially come 
from different excitation conditions in the three sources, and the excitation-abundance 
structure degeneracy can only be strictly broken by observation of additional lines. The 
simplest scenario for explaining our detection is for HNCO to have both an extended and 
a compact origin, however, and this is also supported by the reported excitation charac¬ 
teristics and emission profile of HNCO in other sources (Bisschop et al. 2007). 

In Fig. 5 the signal-to-noise ratio is lower, but it is still clear that CH 3 CHO toward 
NGC7538 IRS9 only has extended emission since none of the IRAM 30m line flux is re¬ 
covered in the SMA spectra. No CH 3 CHO lines are detected in the other two MYSOs in 
the spectral range where IRAM 30m and SMA observations overlap. CH 3 OCH 3 is de¬ 
tected toward NGC7538 IRS9 and AFGL490, and in both cases tentative SMA detections 
suggest that the emission originates in the source centers. Based on the different emission 
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Table 8 : HNCO, CH 3 CHO, and CH 3 OCH 3 lines data from the IRAM 30m spectra and 
HNCO line data from the SMA 2" radius compact region. 


Species 

Freq 

(GHz) 

E up 

(K) 

logA 


Transition 

NGC7538IRS9 

J FdV FHWM 

(Jykms -1 ) (kms -1 ) 

W3IRS5 

$ FdV FHWM 

(Jy km s -1 ) (km s -1 ) 

AFGL490 

$ FdV FHWM 

(Jykms -1 ) (kms -1 ) 

HNCO 

240.876 

112.6 

-3.72 

23 

Hi ,11 - 10 i,io 

2.3+0.6 

13+2 

2 . 1 + 0 .5 

5.5+ 0.5 

< 1.0 

- 

IRAM 

241.704 

239.9 

-3.74 

23 

112,10 -102,9 

< 1.1 

- 

< 1.1 

- 

< 1.0 

- 


241.708 

239.9 

-3.74 

23 

112,9- 102,8 

< 1.1 

- 

< 1.1 

- 

< 1.0 

- 


241.774 

69.6 

-3.71 

23 

llo,n - 10 o,io 

5.9+0.4 

2.3+ 0.3 

4.4+0.7 

3.8+ 0.2 

2.1+0.3 

1.7+ 0.4 


242.640 

113.1 

-3.71 

23 

Hi, 10 - 10 l ,9 

< 1.1 

- 

2.5+0.6 

8.9+ 0.8 

< 1.0 

- 

HNCO 

240.876 

112.6 

-3.72 

23 

Hi ,11 - 10i,io 

< 1.1 

- 

< 1.1 

- 

< 1.0 

- 

SMA 

241.774 

69.6 

-3.71 

23 

llo,n - 10 o,io 

3.2+0.7 

5.8+ 0.6 

4.4+0.9 

5.2+ 0.4 

< 1.0 

- 


242.640 

113.1 

-3.71 

23 

Hi, 10 - 10 l ,9 

1 . 8 + 0.6 

10 + 2 

< 1.1 

- 

< 1.0 

- 

CH3OCH3 

225.599 

69.8 

-3.88 

450 

12 l ,12 - llo.u 

1.7+0.4 

4.2+ 0.4 

0.6+0.3 

2 . 6 + 0.6 

1.5+0.5 

8+2 


240.985 

26.3 

-3.99 

154 

53,3 - 42,2 

<1.1 

- 

<0.8 

- 

<1.3 

- 


241.529 

26.3 

-3.99 

198 

53,2 - 42,3 

<1.1 

- 

<0.8 

- 

<1.1 

- 


241.947 

81.1 

-3.78 

378 

13l,13 - 12o,12 

1.1+0.4 

3.5+ 0.8 

0.4+0.1 

1.2+ 0.3 

1.3+0.5 

6+2 

CH 3 CHO 

223.650 

72.3 

-3.41 

50 

12142-111,!! E 

1.1+0.3 

2.8+ 0.3 

<0.9 

- 

<1.0 

- 


223.660 

72.2 

-3.41 

50 

12i,i2 - lli, 11 A 

1.4+0.3 

3.2+ 0.3 

<0.9 

- 

<1.0 

- 


226.552 

71.4 

-3.39 

50 

12 o,i 2 - H 041 E 

1.5+0.3 

2.9+ 0.3 

<0.9 

- 

<1.0 

- 


226.593 

71.3 

-3.39 

50 

12 o,i 2 - H 041 A 

1.7+0.4 

3.5+ 0.3 

<0.9 

- 

<1.0 

- 


229.775 

61.5 

-4.29 

46 

Hi ,11 - IO 040 A 

<1.0 

- 

<0.9 

- 

<1.0 

- 


230.302 

81.0 

-3.38 

50 

122,11 - H 240 A 

1.0+0.3 

2.5+ 0.3 

<0.9 

- 

<1.0 

- 


230.316 

81.1 

-3.38 

50 

122,11 - H 240 E 

1.0+0.2 

2.3+ 0.3 

<0.9 

- 

<1.0 

- 


242.106 

83.9 

-3.30 

54 

13i,i3 “ 12i,i2 E 

1.1+0.4 

2.6+ 0.5 

<0.9 

- 

<1.0 

- 


242.118 

83.8 

-3.30 

54 

13i,i3'12i,i2 A 

1.4+0.3 

3.7+ 0.5 

<0.9 

- 

<1.0 

- 


244.789 

83.1 

-3.29 

54 

13o,13 “ 12 o,i 2 E 

1.2+0.3 

3.7+ 0.5 

<0.9 

- 

<1.0 

- 


244.832 

83.1 

-3.29 

54 

13o,i3 ‘ 12 o,i 2 A 

1.2+0.5 

4.2+ 0.9 

<0.9 

- 

<1.0 

- 


244.854 

72.3 

-4.19 

50 

12 i,i 2 - llo, 11 E 

<1.0 

- 

<0.9 

- 

<1.0 

- 


patterns, the molecules found in these spectra are classified as follows: CH 3 CCH and 
CH 3 CHO are envelope organics, CH 3 CN and CH 3 OCH 3 are core organics, and CH 3 OH 
and HNCO are intermediate cases with significant core and envelope contributions. 


3.3. Rotational temperatures, and column densities 

The core and envelope classifications based on spatial emission patterns should be re¬ 
flected in the rotational temperatures of the different molecules. Figure 6 shows the rota¬ 
tional diagrams for molecules with enough line detections, i.e., CH 3 OH (extracted from 
IRAM 30m and SMA spectra) and CH 3 CN, CH 3 CCH, following the method described in 
Goldsmith & Langer (1999). The line fluxes from the IRAM observations were converted 
into main beam temperature using the flux-to-antenna temperature conversion factor and 
the beam and forward efficiencies listed online for the EMIR receiver 9 and linearly extrap¬ 
olated for each line frequency (T m t,(K) ^ 0.2 x Flux(Jy)). The line fluxes from the SMA ob¬ 
servation were converted into temperature using the Rayleigh-Jeans approximation with 
a circular beam of 2 /; radius coming from the mask dimension (T m b(K) ~ 1.3 x Flux(Jy)). 
Optically thin emission was assumed based on the low line intensities and lack of asym¬ 
metry in the line profiles. This assumption was verified by the shape of the rotational 
diagram; i.e., flattening or large scatter was observed for lower energy transitions (cf., 
Bisschop et al. 2007). Considering the possibilities of subthermal excitation in the en¬ 
velope, the rotational temperatures are not expected to be the gas kinetics temperatures 
outside of the core. A 10% uncertainty was added to the line-integrated area and is listed 

9 http://www.iram.es/IRAMES/mainWiki/Iram30mEfficiencies 
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239.10 239.15 239.20 239.25 

Frequency / GHz 


Fig. 3: Spectral window with several CH 3 CN and CH 3 CCH lines from the single-dish 
(black lines, 0.2 MHz spectral resolution) and the 2" interferometric data (red line, 
0.8 MHz spectral resolution). 

in the tables to account for the line shapes sometimes deviating from the Gaussian shape 
assumed for the fit. We used the 'linfit' IDL routine to derive the rotational temperatures, 
as well as the column densities, and the routine returned the corresponding uncertainties. 

Table 9 presents the column densities and rotational temperatures derived for these 
molecules using the rotational diagrams in Fig. 6 . The beam-averaged CH 3 OH column 
densities and rotational temperatures derived from the IRAM 30m spectra agree with 
those found by van der Tak et al. (2000), based on JCMT single-dish telescope at higher 

Article number, page 12 of 28 












































Edith C. Fayolle et al.: New constraints on the origins of organics 



241.70 241.75 241.80 241.85 241.90 

Frequency / GHz 


Fig. 4: Spectral window with several CH 3 OH lines from the single-dish (black lines, 
0.2 MHz spectral resolution) and the 2" interferometric data (red line, 0.8 MHz spectral 
resolution). CH 3 OH lines with upper level energies higher than 70 K are marked with a 
star in the NGC 7538 IRS9 to emphasize the increase in SMA/IRAM overlapping for lines 
with higher upper energy levels. 


frequencies. The rotational temperature and column densities derived for CH 3 OH from 
the SMA data are always higher than those derived by the IRAM 30m, which is con¬ 
sistent with the SMA observations probing material closer to the MYSO centers. The 
E-/A- CH 3 OH ratio is consistent with unity within the uncertainties, which agrees with 
Wirstrom et al. (2011). The derived column densities for CH 3 CN from the IRAM data as- 
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242.09 242.11 242.13 241.92 241.94 241.96 

Frequency / GHz Frequency / GHz 


Fig. 5: Spectra of two CH 3 CHO lines at 242.106 GHz and 242.118 GHz (left panel) and 
the CH 3 OCH 3 line at 241.946 GHz (right panel) from the single-dish (black line, 0.2 MHz 
spectral resolution) and 2" interferometric data (red line, 0.8 MHz spectral resolution). 
The three line frequencies are marked by black dotted line. 


sume that the emission is only coming from the 2" radii encompassed by the SMA beam: 
i.e., we apply a dilution factor of 0.16 to account for the SMA extraction mask area (2" ra¬ 
dius) to IRAM beam (5" radius) ratio. This assumption is justified by the CH 3 CN hot-core 
like rotational temperatures of 80-110 K toward the different MYSOs, which are also con¬ 
sistent with the CH 3 OH excitation temperature derived from the SMA spectra. It is also 
consistent with the observed overlap between the IRAM and SMA line fluxes (see Figure 
3). The rotational temperatures of ~50 K obtained for the CH 3 CCH 14-13 ladder from 
the IRAM spectra are consistent with an envelope origin, but suggests that it is mainly 
present in the luke-warm envelope regions rather than in the outermost cold envelope. 

The rotational diagram method assumes that all data can be described by a single 
excitation temperature. To test this assumption for our data, a two-temperature fit was 
explored for the case of methanol. Two-temperature fits of CH 3 OH lines was investigated 
by van der Tak et al. (2000), Leurini et al. (2007), and Isokoski et al. (2013), among others. 
The full results are presented in appendix A, but briefly: both a cold and warm compo- 
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Fig. 6 : Rotational diagrams of CH 3 OH from the single dish data (first column), for CH 3 OH 
from the SMA spectra extracted with a 2" mask (second column), and CH 3 CN from 
the single-dish data (third column), and CH 3 CCH (fourth column) and for the sources 
NGC7538 IRS9 (first row), W3 IRS5 (second row), and AFGL490 (third row). 


nent are recovered from the IRAM 30m data. The derived column densities of the cold 
components are consistent with the single-component fits (within uncertainties), while 
the warm component is consistent with the single-component fit to the SMA data. The 
fit to the SMA line data was not improved by adding a second component, verifying our 
hypothesis that the 2 ,/ mask emission is dominated by a hot component for all sources. 

For FINCO, CFI 3 CFIO, and CFI 3 OCH 3 , no rotational diagrams could be built owing to 
the very small upper-level energy range of the observed transitions, and column densities 
were calculated using the envelope temperature (the CFI 3 OFI IRAM 30m rotational tem¬ 
perature) if the molecule was classified as an envelope molecule and the core temperature 
(the CFI 3 CN rotational temperature) if the molecule was classified as a core molecule, 
and both rotational excitation temperatures if the molecule was classified as intermedi¬ 
ate, i.e. HNCO. As seen in Table 10, the calculated HNCO abundance with respect to 
CFI 3 OH is almost identical regardless of the assumed spatial origin of the line emission. 
For core molecules, the same dilution factor as for CFI 3 CN was applied (see Table 10). For 
molecules with multiple line detections, the column densities were derived by averaging 
the individual column densities found for each detected line and taking the square root of 
the sum of the individual uncertainties squared as uncertainty. Only the IRAM data were 
used to calculate these column densities since these data present a higher signal-to-noise 
ratio. 
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Table 9: Rotational temperatures and column densities for CHsOH, CH 3 CN, and 
CH 3 CCH derived from the rotational diagrams presented in Fig. 6 . B 0 |, s refers to the 
beam radius of the telescope used to obtain the molecular lines (5 "for IRAM and 2"for 
SMA) and B em i ss to the beam size where the line emission is assumed to be coming from. 


Species 

^obs / B em i ss 

NGC7538IRS9 

T ro t (K) N (cm" 2 ) 

Trot (K) 

W3 IRS5 

N (cm -2 ) 

Trot (K) 

AFGL490 

N (cm -2 ) 

CH 3 OH (IRAM) 

5 " / 5 " 

25 ± 2 

9 + 1 x 10 i4 

63+6 

3.2 ± 0.4 x 10 14 

25 + 2 

2.4+ 0.4 x 10 14 

CH 3 OH (SMA) 

2 "/ 2" 

81+ 16 

2.5 + 0.4 x 10 15 

166+ 85 

2.4 ± 0.5 x 10 15 

102 ± 67 

1.2 ± 0.5 x 10 15 

CH 3 CN 

5" IT!' 

111 + 20 

7+ 2x 10 13 

92+23 

2.4 ± 0.8 x 10 13 

164+ 78 

3.6 ± 1.2 x 10 13 

CH 3 CCH 

5" / 5" 

47+ 5 

1.2+ 0.3 x 10 15 

58+8 

7+ 2 x 10 14 

41+ 7 

4.2 ± 1.7 x 10 14 


Table 10: Column densities for HNCO, CH 3 CHO, CH 3 OCH 3 using excitation tempera¬ 
tures from Table 7. B () |, s refers to the beam size of the telescope used to obtain the molec¬ 
ular lines (5 "for IRAM and 2"for SMA) and B em j ss to the beam size from where the line 
emission is assumed to be coming from. 


Species 

®obs / B em i ss 

NGC7538 IRS9 

N (cm~ 2 ) 

W3 IRS5 

N (cm -2 ) 

AFGL490 

N (cm -2 ) 

HNCOext 

5"/ 5" 

4.0 + 1.3 x 10 13 

1.8 ± 0.3 x 10 13 

6.0 + 2.2 x lCF 

HNCOcomp 

5"/ 2" 

1.1 + 0.3 x 10 14 

1.4 + 0.6 x 10 14 

< 4 x 10 13 

CH 30 CH 3 

5"/ 2" 

3.3 + 1.0 x 10 14 

9.7 + 4 x 10 13 

4.5 + 1.9 x 10 14 

CH 3 CHO 

5"/ 5" 

3.1 + 0.4 x 10 13 

< 1.1 x 10 13 

< 1.3 x 10 13 

Table 11: CH 3 CN, CH 3 CCH, HNCO, CH 3 CHO, CH 3 OCH 3 abundances with respect to 
CH 3 OH for their respective spatial origin. 

Species ratio 
/ CH 3 OH 

spatial origin 

NGC7538 IRS9 

W3 IRS5 

AFGL490 

CH 3 CN 

compact 

2.8 + 0.9 x 10 ~ 2 

1.0 ± 0.4 x 10 - 2 

3.0 ± 1.6 x 10 ~ 2 

CH 3 CCH 

extended 

1.3 + 0.4 

2.2 + 0.7 

1.8 + 0.8 

HNCO 

extended 

4.4 + 1.6 x 10 ~ 2 

5.6 + 1.2 x 10 ~ 2 

2.5 + 1.1 x 10 ~ 2 

HNCO 

compact 

4.4 + 1.4 x 10 - 2 

5.8 + 2.7 x 10 ~ 2 

< 3.4 x 10 ~ 2 

CH 3 OCH 3 

compact 

1.3 ±0.5x 10 - 1 

4.0 + 1.8 x 10 - 2 

3.8 + 2.3 x 10 - 2 

CH 3 CHO 

extended 

3.4 +0.6 x 10 - 2 

< 3.5 x IQ ' 2 

< 5.5 x 10 ~ 2 


3.4. Organics in hot cores vs organic-poor MYSOs 

By definition, the organic-poor MYSOs reported in this study have less intense emission 
of complex organic molecules than do line-rich hot cores. The question for this section is 
whether the chemical composition with respect to CH 3 OH is different between the two 
source families. The CH 3 CN, CH 3 CCH, HNCO, CH 3 CHO, and CH 3 OCH 3 abundances 
with respect to CH 3 OH obtained here and presented in figure 11 for the three MYSOs are 
compared to the hot core abundances derived by Bisschop et al. (2007) in Figure 8 . More¬ 
over, three high-mass protostellar objects from Isokoski et al. (2013) have been added. 
These three objects were inferred to have large equatorial structures, but Isokoski et al. 
(2013) found no strong chemical differences in their chemistry compared to the Bisschop 
et al. (2007) hot-core sources. For the three MYSOs sources, the molecular abundances 
with respect to CH 3 OH are calculated using the CH 3 OH column densities derived for 
the envelope if the molecule has been classified as 'envelope' molecules and using the 
CH 3 OH column density derived for the core (SMA-based) in the case of a core molecule. 
For the hot core sources, Bisschop et al. (2007) applied a dilution factor corresponding 
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N[CH 3 OCH 3 ]/N[CH 3 OH]eomp 



N[CH 3 CN]/N[CH 3 0H]comp 



Fig. 7: Gas abundance correlation between organics including upper limits. The black 
crosses are the abundances derived for the MYSOs, the red squares are derived for hot- 
core sources by Bisschop et al. (2007), and the blue diamonds are results for hot-core 
sources from Isokoski et al. (2013). An arbitrary error of 20% has been taken when not 
reported in the two latest studies. 


to the region where T > 100 K for CH 3 OH, CH 3 CN, HNCO, and CH 3 OCH 3 , but not for 
CH 3 CCH and CH 3 CHO. To calculate CH 3 CCH and CH 3 CHO abundances with respect 
to CH^OH, we removed the dilution factor for CH 3 OH applied by Bisschop et al. (2007). 
All other abundances were taken directly from Bisschop et al. (2007). Isokoski et al. (2013) 
could identify a cold and hot methanol emission using single-dish data and the derived 
column densities for both components have been used in figure 7 in the same way as for 
the sources analyzed in the present study. 

The histograms in Figure 8 show that the CH 3 CN, CFI 3 OCFI 3 , and FINCO core abun¬ 
dances with respect to CH 3 OH are similar for the organic-poor MYSOs and the hot-core 
sources. In contrast, the organic-poor MYSOs show higher complex organic envelope 
abundances, i.e., CH 3 CFIO and CFI 3 CCH, with respect to CH 3 OH compared to the hot 
core sources. This difference is most likely due to to our not being able to separate CFI 3 OH 
core and envelope emission in the study by Bisschop et al. (2007), resulting in artificially 
low envelope ratios with respect to CH 3 OH when all CFI 3 OH is implicitly assumed to 
originate in the envelope; in reality, the high excitation temperature of CFI 3 OH in the hot 
core sources suggests that most of it really comes from the core. 

A similar apparent separation between hot core and organic-poor MYSOs are visible 
in log-log correlations of molecular abundances with respect to CFI 3 OFI shown in Fig¬ 
ure 7. Furthermore, there is a clear correlation between envelope molecules CH 3 CFIO 
and CFI 3 CCFI, but this may simply be due to the different abundance derivations of cold 
molecules with respect to CFI 3 OFI for the hot cores and the weak-line MYSOs, rather than 
signifying a chemical relationship. More interestingly, these log-log abundance ratio plots 
show that there is no correlation between the two N-bearing organics CFI 3 CN and FINCO 
over an order of magnitude range. There is also no correlation between the two O-bearing 
complex species CFI 3 OCFI 3 and CFI 3 CFIO, which is consistent with their inferred differ¬ 
ent origins in the organic-poor MYSOs. 
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Fig. 8 : Number of sources versus the logarithm of their gas phase organic ratio over 
methanol with respect to the mean for O-bearing species. The solid filled histograms 
correspond to sources observed and analyzed here, the unfilled histograms correspond to 
sources from Bisschop et al. (2007) and Isokoski et al. (2013). The left panel presents the 
O-bearing species data, while the right panel focuses on the N-bearing molecules. HNCO 
abundances were derived assuming either hot compact emission or cold extended emis¬ 
sion, since its origin does not seem to be consistent between sources. 


3.5. An ice-gas connection? 

The CH 3 OH ice content may be an important factor in whether a hot core chemistry devel¬ 
oped, so we compare the CH 3 OH core column density toward our organic-poor MYSOs 
and hot cores with their CH 3 OH ice abundance (with respect to H 2 O) (see Table 1, Gibb 
et al. 2004). In the hot cores, the majority of the CH 3 OH gas originates in the core and we 
use the derived column densities from Bisschop et al. (2007), where all CH 3 OH emission is 
assumed to originate in the central region where the temperature is higher than 100 K. To 
ensure a fair comparison we calculated the size of the 'hot core region' toward our sample 
using the relation between luminosity and temperature 1 ?t = iook ~ 2.3 x 10 u (^/L/Lq), 
which was shown by Bisschop et al. (2007) to approximate the 100 K radius well toward 
their source sample. We then assumed that all SMA CH 3 OH line flux originate in these 
regions, based on the derived rotational temperatures, and used an appropriate dilution 
factor when the 100 K area is smaller than the 2 // mask used for spectral extraction. Fig- 
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ure 9 presents the resulting column density of hot CH 3 OH gas versus the initial CH 3 OH 
abundance on the grains. It appears that it is primarily the column density of CH 3 OH 
that is different between the line-rich and line-poor sources. No strong correlation with 
the ice content is observed, but more sources would allow the sample to be divided into 
luminosity and total mass bins, removing scatter due to initial physical conditions and 
size of 100 K region. Still, this plot suggests that initial CH 3 OH ice content alone does not 
determine the richness of the MYSO chemistry when correcting for the source luminosity. 
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Fig. 9: CH 3 OH column density in the inner core (calculated area where T > 100 K) versus 
CH 3 OH ice abundance over H 2 O ice in the envelope. The stars present the line-poor 
sources analyzed in this study, while the empty squares are data from Bisschop et al. 
(2007) for hot-core sources. 


The initial ice composition may also affect the complex organic composition in both 
hot cores and organic-poor MYSOs cores and envelopes. Figure 10 presents correlation 
plots between ratios of the N-bearing organics and CFI 3 OFI in the ice and gas phases. 
The two gas-phase N-bearing organics are FINCO and CH 3 CN, and the two ice species 
are OCN and NH 3 . The ice abundances are listed in Table 1 and have been obtained 
by Gibb et al. (2004). For sake of consistency, only the abundances obtained through the 
analysis technique described and performed by Gibb et al. (2004) were used here though 
detailed analysis of specific ice species, taking ice environment and using multiple vibra¬ 
tional bands into account, have been conduced, e.g., by Taban et al. (2003) in the case of 
NH 3 in W 33A. 
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When combining our new observations with data from the literature, a sample of 
seven MYSOs have both ice and gas observations. As seen in Fig. 9, only a fraction of 
them can be used to correlate specific ice ratios, however, because of multiple ice abun¬ 
dance upper limits for many of the sources. For example, W3 IRS5 is not included in any 
of the plots, because of its CFI 3 OFI, OCN and NFI 3 ice upper limits. This means that 
the current data set can only be used to search for tentative correlations or to note gross 
deviations from expected correlations, and not for a proper statistical correlation analysis. 

The top lefthand panel of Figure 10 shows no conclusive correlation between OCN~ 
ice and HNCO in the gas phase with respect to CH 3 OH. Qualitatively, such a correla¬ 
tion is expected since FINCO and OCN~ are linked through efficient thermal acid-base 
chemistry within the ice (e.g., Demyk et al. 1998; van Broekhuizen et al. 2004; Theule et al. 
2011). The lack of a correlation may therefore simply be due to the difficulty determining 
the OCN~ abundance in the ice, so we also explore the correlation between HNCO gas 
and the better constrained NH 3 ice. NH 3 is likely the major source of nitrogen in the ice 
and may therefore be a proxy for the abundance of N-bearing ices in general. It may also 
affect the HNCO/OCN chemistry directly since it is a strong base. The top righthand 
panel of Figure 10 shows that there is indeed a tentative correlation between gas phase 
abundance of HNCO over CH 3 OH with respect to NH 3 . 

The relation between OCN~ in the ice and CH 3 CN in the gas is explored in the bottom 
lefthand panel of Figure 10. These molecules do not appear to be correlated for the four 
sources presented here, despite both containing a CN functional group. Finally, Rodgers 
& Charnley (2001) predict a correlation between CH 3 CN gas a NH 3 ice, but in this lim¬ 
ited sample we find no correlation between the CH 3 CN/CH 3 OH gas ratio versus the 
NH 3 /CH 3 OH ratio in the ice. 

4. Discussion 

4.1. Organic-poor MYSOs versus hot cores 

Previous observations of complex molecules toward MYSOs have generally focused on 
sources with a bright hot core that is responsible for most of the molecular emission. In 
such cases, either interferometric or single-dish observations are sufficient to determine 
complex organic abundances as long as the radius of the evaporation front close to the cen¬ 
tral protostar is known. Single-dish observations combined with the rotational diagram 
technique can also be used to derive abundances of molecules that are predominantly 
present in the outer envelope, since beam-averaged abundances can then be assumed. The 
real difficulty arises for molecules that are distributed throughout the envelope and core. 
For such molecules, single-dish and interferometric observations need to be combined to 
deduce what fraction of the molecular emission originates in the envelope and in the core, 
and then use these fractions to calculate the chemical composition of the two physically 
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Fig. 10: Ice versus gas abundance correlation for N-bearing species with respect to 
CH 3 OH. The crosses are abundances derived for our organic-poor MYSOs, and the red 
squares are the values derived by Bisschop et al. (2007). An arbitrary error of 20% has 
been assumed for the latter values. For the two top plots, the black crosses represent 
the FINCO over CFI 3 OH abundance derived for the compact component, while the blue 
crosses correspond to the HNCO abundance calculation for an extended component. 


and chemically different regions. Based on this study, this class of molecules seems to 
mainly encompass zeroth-generation ices, i.e. CFI 3 OFI and FINCO, but as our sensitivity 
increases, we expect, based on model results, that many classical hot core molecules will 
present a significant envelope emission profile as well (Oberg et al. 2013). 

Using the IRAM 30m and SMA spectra, we could classify several complex organic 
molecules as belonging to the core, envelope, and both. The two envelope molecules, 
CFI 3 CFIO and CFI 3 CCFI, were similarly classified by Bisschop et al. (2007) based on ex¬ 
citation temperatures alone, suggesting that the envelopes around line-poor MYSOs and 
hot cores are chemically similar. In contrast we find that in the line-poor MYSOs, CFI 3 OFI 
and sometimes FINCO have significant emission contributions from the envelope, while 
Bisschop et al. (2007) find that in hot core sources, they have excitation temperatures above 
100 K and were thus classified as originating exclusively in the core region; these sources 
have probably a similar envelope line flux to the one observed for the line-poor MYSOs, 
but in single-dish studies, this emission contribution is drowned out by the hot cores. 
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Overall, the chemistry in the young MYSOs is remarkably similar to what is observed 
in the hot cores, which suggests that they may be hot core precursors. CH 3 CN, CH 3 CCH, 
CH 3 CHO, HNCO, and CH 3 OCH 3 are observed in both kinds of sources at comparable 
abundances with respect to CH 3 OH. CH 3 CH 2 OH and HCOOCH 3 - two typical hot-core 
molecules - are not seen in the organic-poor MYSOs, but typical abundance of these 
molecules with respect to CH 3 OH from Bisschop et al. (2007) are consistent with non¬ 
detections. The hot-core precursor interpretation is also consistent with the observed 
lack of correlation between CH 3 OH core column density or hot-core activity on the initial 
CH 3 OH ice abundance. 


4.2. The ice-gas connection: Observations vs. theory 
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Fig. 11: Model ice versus gas abundance correlation from the MAGICKAL model (Garrod 
2013). Five initial ice abundances (Ml to M5) are used to run the models and derived 
N-bearing abundances with respect to CFI 3 OFI. The black plus signs present the results at 
20 K, the blue crosses are the model results at 30 K, the red squares represents 50 K, and 
the green diamonds are the results at 100 K. 


Regardless or whether the overall hot-core chemistry depends on the initial ice compo¬ 
sition, we expect that the ice composition will have an effect on the chemical composition 
in both hot cores and line-poor MYSOs. This dependence may look very different for com¬ 
plex molecules that form in the gas phase from evaporated ices compared to products of 
complex ice chemistry. That we do not observe a clear trend between NFI 3 in the ice and 
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CH3CN in the gas suggests that the model of Rodgers & Charnley (2001) is missing impor¬ 
tant complex molecule formation pathways. We have therefore used the state-of the-art 
chemical model MAGICKAL (Garrod 2013) to explore the connection between ice and gas 
phase species further. The model uses a rate-equation/modified rate-equation approach, 
treating the gas phase, ice surface, and bulk ice as coupled, but distinct, chemical phases. 

Garrod (2013) produced generic, single-point hot-core models that treated first a cold 
collapse to 10 7 cm” 3 , followed by a warm-up from 8 to 400 K at fixed density, assum¬ 
ing a typical grain size of 0.1 }im and 10 6 surface binding sites. We reran the warm-up 
phase, adopting the medium warm-up timescale of 2 x 10 5 years to reach 200 K (whose 
results appear best to fit various other observational results), but altering the ice abun¬ 
dances prior to warm-up. The original H2O ice abundance is retained, while the CH3OH, 
CH4, HNCO, and NH3 values are varied to mimic the observed ranges in the combined 
line-poor MYSOs and hot cores sample (see Ml to M5 in Table 12). Even though the 
ice abundance ratios correspond to the sources in this study, it is not possible to directly 
compare observations and simulation since key physical parameters such as densities and 
warm up rates of each specific object is not taken into account. The simulations are instead 
used here to investigate chemical trends. The resulting gas-phase abundances of HNCO, 
CH3CN, and CH3OH are reported at different temperatures during warm-up in Table 13. 

Table 12: Initial ice abundances with respect to water used for the five chemical model 
simulation Ml-5 . The three-phase model assumes a H2O abundance with respect to hy¬ 
drogen of 2.66 x 10 7 n [ [ for the ice surface and 1.46 x 10 4 n | [ for the ice bulk. 


Species ratio 

Ml 

M2 

M3 

M4 

M5 

CH3OH/H9O 

4 

11 

14 

4 

10 

CH4/H2O 

2 

2 

2 

2 

2 

HNCO/H2O 

2 

1 

6 

0.5 

2 

NH3/H2O 

15 

6 

15 

6 

6 


Table 13: Gas abundances with respect to hydrogen at various temperatures derived by 
the Garrod (2013) model for the initial ice abundances presented in Table 12 


Model 

Species 

n(20 K)/n H 

n(30 K)/n H 

n(50 K)/n jj 

n(100 K)/n H 

Ml 

CH3OH 

1.6 x 10“ n 

1.1 

X 

10 32 

4.9 x lO”” 11 

3.2 

X 

10~y 


HNCO 

1.7 x 10~ 12 

1.2 

X 

10“ 12 

5.3 x 10" 10 

1.9 

X 

10~ 10 


CH 3 CN 

6.8 x 10” 15 

8.5 

X 

10~ 14 

2.4 x 10~ 10 

3.1 

X 

10-1° 

M2 

CH3OH 

2.0 x 10“ n 

1.9 

X 

10" 12 

1.1 x 10~ 1U 

OO 

YO 

X 

10~ 9 


HNCO 

1.0 x 10~ 12 

1.2 

X 

10~ 12 

4.5 x HT 10 

1.5 

X 

10" 11 


CH 3 CN 

6.8 x 10" 15 

0° 

kO 

X 

10~ 14 

3.9 x 10~ 10 

2.8 

X 

10~ 9 

M3 

CH3OH 

2.2 x lO -1 ' 1 

2.6 

X 

10" 12 

1.1 x 10~ 1U 

7.3 

X 

Ter 9 


HNCO 

4.4 x 10 -12 

1.5 

X 

10" 12 

7.3 x 10" 10 

2.3 

X 

10~ 9 


CH 3 CN 

6.8 x 10” 15 

1.1 

X 

hr 13 

3.7 x 10~ 10 

OO 

CO 

X 

10- 9 

M4 

CH3OH 

1.6 x 10“ n 

9.6 

X 

10~ 13 

5.9 x 10-' 11 

3.7 

X 

Ter 9 


HNCO 

5.4 x 10~ 13 

1.1 

X 

10~ 12 

3.6 x HT 10 

3.6 

X 

10" 11 


CH 3 CN 

6.8 x 10~ 15 

8.1 

X 

10- 14 

2.4 x 10~ 10 

OO 

bo 

X 

i 0 -i° 

M5 

CH3OH 

1.9 x 10“ n 

1.8 

X 

10 -12 

9.4 x 10-' 11 

6.4 

X 

Ter 9 


HNCO 

1.8 x 10 -12 

1.3 

X 

10“ 12 

4.9 x 10" 10 

2.7 

X 

10-n 


CH 3 CN 

6.8 x 10~ 15 

9.2 

X 

IQ - 14 

3.1 x 10~ 10 

2.5 

X 

10-1° 
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The relationships during protostellar warm-up between gas-phase HNCO and CH 3 CN 
with respect to CH 3 OH, and OCN and NH 3 initial ice abundances with respect to CH 3 OH 
ice are shown in Figure 11 for temperatures between 20 K and 100 K. Since the model 
does not treat ion chemistry in the ice, a full conversion rate of HNCO ice into OCN~ 
has been assumed based on the efficient HNCO to OCN~ conversion derived experimen¬ 
tally by Demyk et al. (1998), van Broekhuizen et al. (2004), Theule et al. (2011), among 
others. The resulting complex molecular gas abundances are regulated by a combina¬ 
tion of temperature and initial ice composition. The sensitivity to ice composition varies 
significantly with temperature however, the gas-phase HNCO/CH 3 OH ratio, for exam¬ 
ple, barely changes with ice composition when the grains are sitting at 30 K, because of 
the very limited sublimation at this temperature. To predict the complex chemistry thus 
clearly requires knowing both the temperature structure and the initial ice composition of 
a source. 

Based on the observational and theoretical results that most HNCO and CH 3 CN emis¬ 
sion start at high temperatures, we focus on the predictions at 100 K. At this temperature, 
the HNCO gas vs. methanol content is, as expected, correlated to the initial amount of 
OCN~ over methanol in the ice. This result suggests that with more sources and/or bet¬ 
ter constraints on OCN~ ice abundances, a clearer correlation in the observed data should 
appear as long as the model captures the dominant HNCO formation/destruction path¬ 
ways. The observed tentative correlation between gas-phase HNCO and NH 3 ice is con¬ 
sistent with model predictions, except for the M3 run. In the M3 model, the high absolute 
abundance of HNCO ice results in a longer HNCO desorption time scale, which shifts the 
abundance peak of HNCO to higher temperatures and results in the high [HNCO]/ [CH 3 OH] 
ratio at 100 K; at higher temperatures, M3 no longer deviates from the trend. To fully ex¬ 
plore the effects of NH 3 and OCN “ on the final complex organics abundances clearly 
requires a much larger grid of models that covers all possible combinations of ice abun¬ 
dances as well as investigating the temperature dependences of the complex chemistry. 

In contrast to what has been proposed by Rodgers & Charnley (2001), the abundance 
of CH 3 CN with respect to CH 3 OH does not correlate with either the NH 3 ice content in 
the MAGICKAL code output or with the cyanide ice-related species OCN~. This agrees 
with the observational results. In Rodgers & Charnley (2001) and (Garrod 2013), CH 3 CN 
forms mainly through radiative association reaction in the gas phase between CH^ and 
desorbing HCN giving CH 3 CNH + . The correlation between CH 3 CN and NH 3 predicted 
by Rodgers & Charnley (2001) comes from a cycled production of HCN from NH 3 . The 
latter is, however, not observed in MAGICKAL, which explains the lack of correlation 
between CH 3 CN and NH 3 . 

In summary, there is some encouraging tentative agreement between model predic¬ 
tions and observations. To directly compare models and observations requires, however. 
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that the appropriate model results are mapped onto the temperature-density profiles of 
individual sources, since both ice composition and temperature are shown to strongly af¬ 
fect the complex chemistry (Oberg et al. 2013). Thus, to draw any general conclusions 
requires a large sample of spatially-resolved gas-phase observations, along with ice ob¬ 
servations of the same object. As shown here, organic-poor high-mass protostars contain 
detectable amounts of complex organic material and present a similar chemistry to bright 
hot cores. Most massive YSOs with existing ice observations could therefore be used to 
expand the sample of sources. 

5. Conclusions 

We detected complex organic molecules CH 3 CN, CH 3 CCH, CH 3 CHO, and CH 3 OCH 3 
together with HNCO and CH 3 OH, toward three massive YSOs without any previous evi¬ 
dence of hot-core chemistry activity. Using a combination of single-dish and interferome¬ 
try observations, we found that CH 3 CN and CH 3 OCH 3 emission originates in the central 
core region, CH 3 CHO and CH 3 CCH in an extended envelope, and CH 3 OH and, some¬ 
times, HNCO have both envelope and core emission components. The inferred molecular 
emission locations are consistent with rotational temperatures derived from the single¬ 
dish observations, except for CH 3 OH, where single-dish data are dominated by the enve¬ 
lope. 

The high-temperature abundances of complex organics with respect to CH 3 OH are 
indistinguishable for the organic-poor MYSOs and the sample of hot core sources from 
Bisschop et al. (2007) and Isokoski et al. (2013). The envelope chemistry also seems similar 
for both kinds of sources, but this analysis is limited by a lack of CH 3 OH envelope data 
toward hot core sources. No strong correlation between initial CH 3 OH ice abundance and 
hot CH 3 OH gas column density close to the central object was observed. 

The NH 3 ice abundances seem to affect the HNCO/CH 3 OH gas-phase abundances. 
This relationship is reproduced by the MAGICKAL astrochemical code, assuming fiducial 
collapse and warm-up rates and initial ice compositions that span the observed range. 

More sources with both ice and gas data are required to settle how ice abundances 
affect complex molecule distributions around MYSOs. This could be achieved by col¬ 
lecting mid-IR ice spectra (using SOFIA for example), and performing spatially resolved 
millimetric observations of sources with no detected hot-core molecules. 
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Appendix A: Two-component fit of the CH3OH rotational diagrams 

The IRAM 30m and SMA CH 3 OH line data have been analyzed further to explore whether 
the data is better fit by two distinct temperature distributions than by the one assumed in 
the rotational diagrams in the main body text. In fact, two Boltzman distributions would 
be expected for the IRAM 30m data, since these spectra include both the envelope and the 
warm core seen in the SMA spectra. Figure A.l presents the two-component rotational 
diagrams for the three MYSOs. The fits were performed using the IDL routine MPFIT 
with the initial temperature guesses of 20 K and 100 K and initial column densities guesses 
of 5 x 10 14 cm 2 and 1 x 10 15 cm 2 . 

For all three MYSOs the data is well fit by two components, but the uncertainties in the 
derived excitation temperatures and column densities are very large, demonstrating that 
the signal-to-noise ratio of this data set is not sufficient for a quantitative two-component 
analysis. Still, it is clear that in each case, the cold component has a slightly lower excita¬ 
tion temperature than the single-component fit. The derived column densities agrees (tak¬ 
ing the large uncertainties into account) between the two fits, supporting our assumption 
that the IRAM 30m CH 3 OH spectra are dominated by the envelope. Each fit also results 
in a warm component. While the excitation temperatures have large uncertainties, they 
are all consistent with those derived from the SMA data, further supporting the conclu¬ 
sions based on that data set. Owing to these large uncertainties for the column densities 
and the hot component temperature, single-temperature fitting was used for the quanti¬ 
tative analysis in the paper. Higher signal-to-noise ratio single-dish data could clearly be 
used, however, to simultaneously constrain temperature and column densities of the two 
components, limiting the need for high-spatial-resolution observations for some sources. 

We also attempted to fit the CH 3 OH SMA spectra with two temperature components 
using the same fitting routine and initial guesses. In each case, the outcome was that a 
single-component fit the data as well as two components; i.e., it was not possible to dis¬ 
tinguish multiple components with different temperatures within the SMA masks. This 
confirms that the SMA beam samples only the hot CH 3 OH component and filters out the 
extended cold emission. 
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